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ABSTRACT 



In this research, the development of a fast numerical method was performed in order to solve the Option Pricing 
problems governed by the Black-Scholes equation using an accelerated genetic algorithm method. 
Where the Black-Scholes equation is a well known partial differential equation in financial mathematics. A discussion of 
the solutions was introduced for the linear Black-Scholes model with the European options (Call and Put) analytically and 
numerically ,with and without transformation to heat equation .Comparisons of the presented approximation solutions of 
these models with the exact solution were achieved. 

KEYWORDS: Linear Black-Scholes Equation, American and European Options, Genetic Algorithm 
1. INTRODUCTION 

Many mathematical models have been proposed to describe the evolution of the value of financial derivatives. 
Financial models were generally formulated in terms of stochastic differential equations, [1]. Linear Black-Scholes 
equation the model was developed by Fischer Black and Myron Scholes [2], as an important equation in the financial 
mathematics. 



Where S:= S(t) > 0 and t e (0, T ), provides both the price for a European option and a hedging portfolio that 
replicates the option assuming that, [4]. 

• The price of the asset price or underlying derivative S(t) follows a Geometric Brownian motion W(t), meaning that 
S satisfies the following stochastic differential equation (SDE). 

dS(t) = [iS(t)dt + aS(t)dW(t). 

• The trend or drift /x (measures the average rate of growth of the asset price), the volatility a (measures the standard 
deviation of the returns) and the riskless interest rate r are constant for 0 <t <T and no dividends are paid in that 
time period. 

Equation (1) can be transformed into the heat equation to find the price of the option analytically, [1]. In this 
research, the application of the accelerated Genetic Algorithm method was taken into account to solve the linear 
Black-Scholes model for European options. Studying (1) for an American options would be redundant, since the value of 
an American options equals the value of a European options if no dividends are paid and the volatility is constant, [6], [7], 



0 = V t + -a 2 V ss + rSV s - rV 



(1) 



[8], [9], [10]. 
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2. THE GENETIC ALGORITHMS 

The principles of genetic algorithm are discussed in previous paper [15]. Where The components of the genetic 
algorithm, [11, 12] are: 

• The initial population. 

• The fitness function. 

• The genetic operators. 

2.1. Evolutionary Algorithm Used for Coding Production 

To produce fitness function we used Backus-Naur grammar form (BNF), [13, 14, 15]. 
Where the grammar into as the following: 

• Reading an element from the chromosome (with value V). 

• Selecting the rule according to the scheme: 
Rule = V mod NR (1) 

Where NR is the number of rules for the specific non-terminal symbol And the grammar shown in Table 1 below. 
If the chromosome g = [7 9 4 14 28 10 12 2 17 15 6 11 10 24 11]. Table 2 show how the function is produced by the 
grammar. Thus the grammatical evaluation can found in, [16, 17, 18]. 



Table 1: The Grammar of the Proposed Method 



S::=<expr> 




<expr> ::= <expr> <op> <expr> 


(0) 


1 ( <expr> ) 


(1) 


1 <func> ( <expr> ) 


(2) 


l<digit> 


(3) 


Ix 


(4) 


iy 


(5) 


Iz 


(6) 


<op> ::= + 


(0) 


1- 


(1) 


1 * 


(2) 


1/ 


(3) 


<func> ::= sin 


(0) 


Icos 


(1) 


lexp 


(2) 


Hog 


(3) 


<digit> ::= 0 


(0) 


1 1 


(1) 


12 


(2) 


13 


(3) 


14 


(4) 


15 


(5) 


16 


(6) 


17 


(7) 


18 


(8) 


19 


(9) 
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Table 2: Illustrate Example of Program Construction 



String 


Chromosome 


Operation 


<expr> 


/ 2 y 4 12 2 20 1 2 j 1 4 


/moa/=0 


<expr><op><expr> 


2 9 4 12220 12 314 


2mod 7=2 


Fun(<expr>)< op><expr> 


94 12220 12 31 4 


9mod4=l 


cos(<expr>)< op><expr> 


412 2 2012314 


4mod7=4 


cos(x)<op><expr> 


12 2 20 12 31 4 


12mod4=0 


cos(x) + <expr> 


2 20 12 31 4 


2mod7=2 


cos(x) + Fun(<expr>) 


20 12 31 4 


20mod4=0 


cos(x) + sin(<expr>) 


12 314 


12mod7=5 


cos(x) + sin(y) 







2.2. Technique of the Proposed Method 

The proposed method has the following phases: 

• Initialization. 

• Fitness evaluation. 

• Genetic operations. 

• Termination control. 

2.2.1. Initialization 

The value of mutation rate and selection rate are stated, [13, 15]. The initialization of every chromosome is 
performed by randomly selecting an integer for every element of the corresponding vector. 

2.2.2. Fitness Evaluation 

Expressing the Partial differential equation in the following form: 

f \ x ' y> Yx y) ' ^ y) ' ( x> y) ' ( X) y) J = 0 ' x G * i] ' y e yi] 

The associated boundary conditions are expressed as: 

"fro, y ) = f 0 (y) ,u(.*i .y ) =fi(y)M* ,y„) = g 0 (y)^(x .y-j = g ± (y) 

The steps for the fitness evaluation of the population are the following: 

• Choose N equidistant points in the box [jEo.XjJ X [y^yj , N x equidistant points on the boundary at x = x 0 and at 
x = X] , N y equidistant points on the boundary at y = y 0 and at y = yt 

• For every chromosome i: 

o Construct the corresponding model M(x ,y), expressed in the grammar described earlier, 
o Calculate the quantity 

f ( Wi ) = Sf= o Cf - y; - ^ w i - y; ) ^ w i ( x j >yj)>^ M i i x ) >yj)>j~i M i t x j ' y; ) ) 2 

o Calculate an associated penalty P, (M t ) . The penalty function P depends on the boundary conditions and it 
has the form: 

www.tiprc.org editor@tjprc.org 



36 



Eman Ali Hussain & Yaseen Merzah Alrajhi 



o Calculate the fitness value of the chromosome as: 



v t = E{M t ) + P ± (M0 + P 2 (M i ) + P 3 (M0 + P 4 (M t ) 



2.2.3. Genetic Operators 

The genetic operators that are applied to the genetic population are the initialization, the crossover and the 
mutation. A random integer of each chromosome was selected to be in the range [0.255]. The parents are selected via 
tournament selection, i.e: 

• First, create a group of K >= 2 randomly selected individuals from the current population. 

• The individual with the best fitness in the group is selected, the others are discarded. 

The final genetic operator used is the mutation, where for every element in a chromosome a random number in the 
range [0, 1] is chosen, [15]. 

2.2.4. Termination Control 

Creating new generation required for application genetic operators to the population in order to find the best 
chromosome having better fitness or whenever the maximum number of generations was obtained. 



3. TECHNICAL OF THE ACCELERATED METHOD 

To make the method is faster to arrived the exact solution of the partial differential equations by the following: 



• Insert the boundary conditions of the problem as a part of chromosomes in the our population of the problem, the 
algorithm gives the exact solution or approximate solution in a few generations. 

• Insert a part of exact solution (or particular solution) as a part of a chromosome in the population, find the 
algorithm that gives an exact solution in a few generations. 

• Insert the vector of exact solution (if exist) as a chromosome in the our population of the problem, the algorithm 
gives the exact solution in the first generation. 



4. EXACT SOLUTION OF THE LINEAR BLACK - SCHOLES EQUATION 



To solved equation (1) in closed form [19, 20]: 




Now 
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2 St 



av _ 




av as 








dt ~ 


dz at 


as at 


dV 




k av 




K^ = 




dS 


as 


S dx 



. , a 2 v „ a 2 v k r a 2 v av, 
And = K — = — \— - — 1 
as 2 as 2 s 2 l ax 2 ax 1 



Substituting these derivatives in (1) one gets 

dv d 2 v (1r .\av 2r 
= [~~ 1 la 2 V 

OT OX L \O l / OX <7 Z 

Letting — = 6 then the above equation can be written as : 



Ji = £f| +(0 _l)Ji_ 017 (2) 

or ax L ax 

Now letting X = -(G - 1) ,v = -(6 + 1) = A + 1 
To get v 2 = X 2 + 0 
And then 

v(x,t) — e~ Xx ~ yl " c u{x,T) , 
Therefore; 

dv _ e -^x-v 2 T^_ v 1^ X T) + ^ e -^-" 2 T u ( X)T ) - e - Xx - v2r u(x,T)[(-V 2 )u + ^1 

OT OT OT 

!z ==e -WT u( )[ _ PLU + ^ 

dx dx 1 

a 2 v -Ax-u 2 Tn2 m 3u , ° 2 "l 

ox 2 L ox dx 2j 

Inserting these into equation (2) and dividing by e~ Xx ~ v2r , obtaining: 



(3) 



The initial and boundary conditions for the European Call and Put options are respectively 

u{x, 0) = (e a+1)x - e Xx ) + as xeR, 

u{x, r) = 0asn -co 

u{x, t) = ( e W+D*+^ _ e ix+i 2 T^ asx^oo 
And 

u(x, 0) = [e Xx - e ^ x+1)x y as xeR, 

u(x, t) = 6 ^ 2 x+Xt as x _ Q0 

u{x, t) = 0 as x -* oo 
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Thus the Black-Scholes equation reduces to the heat equation. Now considering the problem (3) with initial 
condition: 

For Call Option, 

uO,0) = g(x) = max{(e a+1) * - e Ax ),0) , 



And for Put Option 

u(x, 0) = g{x) = max {(e Xx - e a+1)x ), 0} , 

Applying the Fourier transformation with respect to x in equation (3) where it was solved in [22] yields; 
Call Option 

V(S, t) = Kv{x,t) = Ke- Xx - v2r u(x,T) 
Therefore 

V C (S, t) = 50(dl) - Ke- r ^0(d2) 
Where 

dl HIH^)^ 

a^T-t 

And 



dl 



Kf)+K)^) 



dl - a^T - t 



This gives the values for the European Call option where <f>(d) the standard normal distribution of the function is 
d. V(S,t) can be seen in Figure 1. The corresponding parameters are K = 100, a =0.2, r =0 .1, T = 1 year, dt =0 .001. 

Put Option 

A similar calculation for the European Put option and the pay-off function for the Put option in absence of 
dividend is 

V P (S,t) = Ke- r ^0(-d2) -S0(-dl) 

Figure 1 is related to the corresponding parameters K = 100, a =0.2 , r =0 .1, T = 1 year , dt =0 .001. 





(a) 



(b) 



Figure 1: Exact solutions of Eq. (1), (a) Call Option, (b) Put Option 
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5. APPLICATIONS OF ACCELERATED GENETIC ALGORITHM 

To find numerical solutions of Linear Black-Scholes model and the transformed model we applied the accelerated 
genetic algorithm method. The crossover rate was set to 75% (that is, replication rate equal to 25%) and mutation rate was 
set to 1%. Each experiment to find the numerical solution was performed 20 times. The population size was set to 100 and 
the length of each chromosome to 50. The size of the population is a critical parameter. The generations number set to 
100 iterations for each experiment. The function (randi) in Matlab R2010a used to generate the initial population. 

5.1. Numerical Solution by Accelerated Genetic Algorithm of Transformed Equation 

In this section an attempt was made to solve the model numerically. The accelerated genetic algorithm scheme 
was recalled to solve the transform problem; the heat equation. Then backward substitution reduced the solution of 
equation (1). 

Call Option 

Considering the model 
u T = u xx ,x £ R ,0 < t < f 

With the initial and boundary conditions for call option as 
u(x, 0) = (e^ +1)x - e Ax ) + as xeR, 
u(x, t) = 0osi-> —CO 



The approximate solutions was obtained as: 
u(x, t) = e~ T sinx at generation 20 
And 

Gp4(x,T) — (x + xe~ T )/2 at generation 4 , as a trail solutions. 
Then, the approximation solutions of (1) are: 



u(x, t) 



) as x -> oo 




And 



Gp4(S,t) = 0.5S(lng) + mg) 




These solutions and comparisons of them with the exact solution of eq. (1) shown in Figure 2: 



www.tiprc.org 



editor@tjprc.org 



40 



Eman Ali Hussain & Yaseen Merzah Alrajhi 



V(S l 1)=S.sintln(S/hq),e'-ly2 r (T-t)) 





(a) V(s, t) 



(b) Gp4(S, t) 




(c) Exact Solu. of (1) and V(S, t) (d) Exact Solu. of (1) and Gp4 
Figure 2: Trail Solutions and There Comparisons of them with the Exact Solution of Eq. (1) (call Option) 

Put Option 

The related model is; 
u T = u xx ,x £ R ,0 < t < f 

With the initial and boundary conditions for Put option as 
u(x, 0) = [e Xx - e^ A+1)x ) + as xeR, 
u(x, t) = e A2x+Xr as x -> — oo 
u(x, t) = 0 as x -* oo 
And the approximate solutions are; 

u(x,t) = -1 + ^e~ x (e T/8 + 3e T/4 ); at generation 20 , and 
Gp4(x,r) — — 1 + e~ x ~ x ; at generation 4 

The back substitution of the variables transformation can be used to get the solution of a problem related to the 
partial differential equation (1). The numerical solutions are: 



V(S,t) = -S + -(e 



K ^(^C-t))/ 8 _i_ Q„(> 2 (r-t))/4 



; and 



Gp4(S,t) = -S + Ke 
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(c) Exact Solu. of (1) and Gp4 (d) Exact Solu. of (1) and V(S, t) 

Figure 3: Trail Solutions and Comparisons with the Exact Solution of eq. (1); Put Option 

5.2. Numerical Solution of Original Equation by Accelerated Genetic Algorithm 

In this section a trial was achieved to solve the original model numerically. The use of accelerated genetic 
algorithm scheme was recalled to solve the original equation (1). 

0 = V t + \a 2 V ss + rSV s - rV 

For call and put options with boundary and initial conditions 
1/(0, t) = 0/or 0 < t < T 
V(S, t)~S - Ke- r( - T -V as S -> oo 
V(S, T) = (S- K) + for 0 < S < oo 
And 

V(0, t) = Ke r(T ~V forO<t<T 

V(S, t)~OasS^>oo 

V(S, T) = (K - S) + for 0 < S < oo 
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For Call Option 

Gp5(5, t) = V(S, t) = S- 100e"^ (7, " £) 

as a trail solution at generation 5 , as well as 

Gp7(S, t) = V(S, t)=S- 81e-° 02 ^ 
as the trail solution at generation 7. 
For put option the above solution satisfied such that: 

Gp5'(S, t) = V(S, t) = 100e"^ r_t) - 5 
And Gp7\S, t) = V{S, t) = Q\ e -°™(T-t) _ s 
These solutions are represented by in Figure 4. 




(b) Gp5(S, t) (Call Option) (b) Gp7(S, t) (Call Option) 




s 



(c) Gp5'(S, t) Put Option (d) Gp7'(S, t) Put Option 

Figure 4: Numerical Solutions of Original Equation without Transformation 

6. COMPARISON THE STUDY 

Comparison all solutions for European Options (call and put) corresponding parameters K = 100, a =0. 2 , r =0 
.1, T = 1 year , S=0-150 for both call and put options: 

6.1. Call Option 

In Table 3 and Table 4 below we show that the errors between these solutions, and the difference between these 
errors are shown Figure 5-a. 



Impact Factor (JCC): 4.2949 



Index Copernicus Value (ICV): 3.0 



Accelerated Genetic Algorithm Solution of Linear Black - Scholes Equation 



43 



Table 3: Comparison of Solutions for European Call Option 
(Exact Solution with Numerical Solutions of Transformed Eq) 





t 


Exact Solu. 
ofEq.(l) 


V(S, t) of 
Transformed 

Fn <"V1 


Gp4(S, t) for 
Transformed 


Error = 
Abs(l^xact- 


Error = 
Abs(Exact-Gp4) 


a 
U 


1 
1 


A 
U 


A 
U 


a 
u 


A 

u 


A 
U 


1 GO/IT 


U.V4/3 


A 
U 


A 
U 


A 
U 


A 
U 


A 
U 


1 C TOO/I 

1 j. /oV4 


A QftAH 

u.oy4/ 


A 
U 


A 
U 


A 
U 


U 


U 


Zj.oo4Z 


U.84Z1 


o nno Q 1 a 

o.Ul /oe-14 


U 


U 


o mc 1 a 
o.Uzh-14 


o mc i a 
o.UZb-14 


1 1 C7QO 

M.j toy 


U. /o94 


1 1 A A C„ AA 

1. 144je-UV 


0 


n 
U 


1 1 /IE AA 


1 1 /IE AA 


Q A AH1£. 


U. /JOO 


C AA/CA,-. AT 


0 


U 


A r- AT 

ob-U/ 


ofc-U/ 


A H A 


U.oo4Z 




0 


n 
U 


/I O/IC AC 

4.94n-Uj 


/I C\AC AC 

4.y4ri-UD 


CC 16.1. 1 


n 1 c 


A A A 1 O 

U.UUlz 


A 
U 


A 
U 


A AA 1 1 

U.UUlz 


A AA 1 T 

U.UUlz 


/TO 1 CHQ 

63.1 J /o 


U.J /oV4 


A A1 C 1 

U.U1M 


A 
U 


A 
U 


A A1 C 1 

U.U1 jl 


A A1 C 1 

U.U1 jl 


/ I.Ujzo 


a. 


A 1 AOA 

U.lUzV 


A 
U 


A 
U 


A 1 ATA 
U.1U29 


A 1 AOA 

u. iuzy 


78 QzL7^ 


O zL7^ 




A 


A 
\J 






86.8421 


0.4210 


1.5177 


0 


0 


1.5177 


1.5177 


94.7368 


0.3684 


3.8485 


0 


0 


3.8485 


3.8485 


102.6315 


0.3157 


7.9136 


2.6488 


2.6575 


5.2648 


5.2561 


110.5263 


0.2631 


13.6821 


10.9854 


11.0328 


2.6967 


2.6493 


118.4210 


0.2105 


20.5733 


19.8432 


19.9801 


0.7301 


0.5932 


126.3157 


0.1578 


27.8843 


29.1493 


29.4627 


1.265 


1.5784 


134.2105 


0.1053 


35.2576 


38.8408 


39.4485 


3.5832 


4.1909 


142.1052 


0.0526 


42.6302 


48.8626 


49.9092 


6.2324 


7.279 


150 


0 


50 


59.1669 


60.8197 


9.1669 


10.8197 



Table 4: Comparison of Solutions for European Call Option 
(Exact Solution with Numerical Solutions of Eq. 1) 



s 


t 


Exact Solu. of 
Eq.(l) 


Gp5(S,t) of 
Eq.(l) by Our 
Method 


Gp7(S,t) of 
Eq.(l) by Our 
Method 


Error = 
Abs(Exact- 
Gp5) 


Error = 
Abs(Exact- 
Gp7) 


0 


1 


0 


0 


0 


0 


0 


7.8947 


0.9473 


0 


0 


0 


0 


0 


15.7894 


0.8947 


0 


0 


0 


0 


0 


23.6842 


0.8421 


8.0178e-14 


0 


0 


8.02E-14 


8.02E-14 


31.5789 


0.7894 


1.1445e-09 


0 


0 


1.14E-09 


1.14E-09 


39.4736 


0.7368 


5.9960e-07 


0 


0 


6E-07 


6E-07 


47.3684 


0.6842 


4.9382e-05 


0 


0 


4.94E-05 


4.94E-05 


55.2631 


0.6315 


0.0012 


0 


0 


0.0012 


0.0012 


63.1578 


0.57894 


0.0151 


0 


0 


0.0151 


0.0151 


71.0526 


0.5263 


0.1029 


0 


0 


0.1029 


0.1029 


78.9473 


0.4736 


0.4636 


0 


5.2687 


0.4636 


4.8051 


86.8421 


0.4210 


1.5177 


0 


12.3838 


1.5177 


10.8661 


94.7368 


0.3684 


3.8485 


0 


19.4906 


3.8485 


15.6421 


102.6315 


0.3157 


7.9136 


2.9468 


26.5891 


4.9668 


18.6755 


110.5263 


0.2631 


13.6821 


10.7891 


33.6792 


2.893 


19.9971 


118.4210 


0.2105 


20.5733 


18.6313 


40.7607 


1.942 


20.1874 


126.3157 


0.1578 


27.8843 


26.4735 


47.8337 


1.4108 


19.9494 


134.2105 


0.1053 


35.2576 


34.3157 


54.8979 


0.9419 


19.6403 


142.1052 


0.0526 


42.6302 


42.1578 


61.9534 


0.4724 


19.3232 


150 


0 


50 


50 


69 


0 


19 



6.2. Put Option 

In Table 5 and Table 6 below we show that the errors between these solutions, and the difference between these 
errors are shown Figure 5-b. 
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Table 5: Comparison of Solutions for European Put Option 
(Exact Solution with Numerical Solutions of Transformed Eq) 







Exact Solu. 


V(S, t) of 
Transformed 


Gp4(S ,t) for 
Transformed 


Error = 
Abs(Exact- 


Error= 
Abs(Exact- 






90.4837 


100.4385 


98.0198 


9.9548 


7.5361 




U.V4/3 


83.0664 


92.5206 


90.2283 


9.4542 


7.1619 


1 ^ TOO/1 


U.oV4/ 


75.6517 


84.6027 


82.4369 


8.951 


6.7852 


Zj.0o4z 


O C/IO 1 


68.2395 


76.6849 


74.6456 


8.4454 


6.4061 


31. j /oV 


n HQCiA 


60.8299 


68.7670 


66.8545 


7.9371 


6.0246 




U. /3oo 


53.4228 


60.8492 


59.0634 


7.4264 


5.6406 


zL7 ^6RzL 


U.Oo^fZ 


46.0183 


52.9313 


51.2724 


6.913 


5.2541 


55.2631 


0.6315 


38.6176 


45.0135 


43.4816 


6.3959 


4.864 


63.1578 


0.5789 


31.2321 


i '—i r\r\ n 

37.0957 


o t~\f\ ct 

35.6908 


5.8636 


4.4587 


71.0526 


0.5263 


23.9232 


29.1779 


27.9002 


5.2547 


3.977 


78.9473 


0.4736 




zl.zoUU 




4.3 /Ul 


i.Zlyo 


86.8421 


0.4210 


10.5525 


13.3422 


12.3193 


2.7897 


1.7668 


94.7368 


0.3684 


5.4945 


5.4244 


4.5290 


0.0701 


0.9655 


102.6315 


0.3157 


2.1735 


0 


0 


2.1735 


2.1735 


110.5263 


0.2631 


0.5585 


0 


0 


0.5585 


L 0.5585 


118.4210 


0.2105 


0.0689 


0 


0 


0.0689 


0.0689 


126.3157 


0.15789 


0.0020 


0 


0 


0.002 


0.002 


134.2105 


0.1052 


1.9454e-06 


0 


0 


1.95E-06 


1.95E-06 


142.1052 


0.0526 


0 


0 


0 


0 


0 


150 


0 


0 


0 


0 


0 


0 



Table 6: Comparison of Solutions for European Put Option 
(Exact Solution with Numerical Solutions of Eq. (1)) 





t 


Exact Solu. 
ofEq.(l) 


Gp5'(S, t) of 
Eq.(l) by 

Our 
Method 


Gp7'(S, t) of 
Eq.(l) by Our 
Method 


Error = 
Abs(Exact- 
Gp5') 


Error = 
Abs(Exact- 
Gp7') 


0 


1 


90.4837 


99.0049 


79.396 


8.5212 


11.0877 


7.8947 


0.9473 


83.0664 


91.1623 


71.5849 j 


8.0959 


11.4815 


15.7894 


0.8947 


75.6517 


83.3197 


63.7739 


7.668 


11.8778 


23.6842 


0.8421 


68.2395 


75.4772 


55.963 


7.2377 


12.2765 


31.5789 


0.7894 


60.8299 


67.6346 


48.1521 


6.8047 


12.6778 


39.4736 


0.7368 


53.4228 


59.7921 


40.3413 


6.3693 


13.0815 


47.3684 


0.6842 


46.0183 


51.9497 


32.5307 


5.9314 


13.4876 


55.2631 


0.6315 


38.6176 


44.1072 


24.7201 


5.4896 


13.8975 


63.1578 


0.5789 


31.2321 


36.2648 


16.9096 


5.0327 


14.3225 


71.0526 


0.5263 


23.9232 


28.4224 


9.0992 


4.4992 


14.824 


78.9473 


0.4736 


16.8899 


20.58 


1.2888 


3.6901 


15.6011 


86.8421 


0.4210 


10.5525 


12.737 


0 


2.1845 


10.5525 


94.7368 


0.3684 


5.4945 


4.8954 


0 


0.5991 


5.4945 


102.6315 


0.3157 


2.1735 


0 


0 


2.1735 


2.1735 


110.5263 


0.2631 


0.5585 


0 


0 


0.5585 


0.5585 


118.4210 


0.2105 


0.0689 


0 


0 


0.0689 


0.0689 


126.3157 


0.15789 


0.0020 


0 


0 


0.002 


0.002 


134.2105 


0.1052 


1.9454e-06 


0 


0 


1.95E-06 


1.95E-06 


142.1052 


0.0526 


0 


0 


0 


0 


0 


150 


0 


0 


0 


0 


0 


0 
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1 3 5 7 3 11 13 15 17 19 

No. of iterations No. of iterations 

(a) Call Option (b) Put Option 

Figure 5: Difference Errors (a) Call Option, (b) Put Option 

7. CONCLUSIONS 

In this study the accelerated genetic algorithm method was used to obtaining a wide class of approximations 
solutions for the linear Black-Scholes model with transformed to heat equation and without transformed, and then choosing 
the best numerical solution among these solutins. One can noted that this method has general utility for applications. 
Finally, the accelerated Genetic Algorithm has convenient used for approximating the analytical solutions for this model. 
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